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We propose in this paper a quintom model of dark energy with a single scalar field cf> given by the 
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lagrangian £ = -V{4>)^/1 - a'V„<j)\7^(j> + (i'<f>ncf). In the limit of f3' ^0 our model reduces to the 
^-j. effective low energy lagrangian of tachyon considered in the literature. We study the cosmological 

evolution of this model, and show explicitly the behaviors of the equation of state crossing the 
| cosmological constant boundary. 

m ■ 

(N ; I. INTRODUCTION 

CO ■ 

The current data from type la supernovae, cosmic microwave background (CMB) radiation, and other cosmological 
\Q ' observations [J d, d, [J] have provided strong evidences for a spatially flat and accelerated expanding universe at the 
present time. Within the framework of the standard cosmology, this acceleration can be understood by introducing 
■ a mysterious component, dubbed dark energy (DE). The simplest candidate for DE is a minor positive cosmological 
constant, but it suffers from the fine-tuning and coincidence problems. As a possible solution to these problems 
various dynamical models of DE have been proposed, such as quintessence. In the recent years with the accumulated 
astronomical observational data it becomes possible to probe the current and even early behavior of DE. Although 
the current fits to the data in combination of the 3-year WMAP[f|, the recently released 182 SNIa Gold sampleQ 
and also other cosmological observational data show the consistence of the cosmological constant, it is worth noting 
that the dynamical DE models are not excluded and a class of dynamical models with equation of state across — 1, 
dubbed quintom, is mildly favored (for recent references see e.g. 0, |p). 

Theoretically it is a big challenge to the model building of the quintom dark energy. The Ref. @ is the first paper 
pointing out this challenge and showing explicitly the difficulty of realizing w crossing over — 1 in the quintessence and 
. ' phantom like models. In general with a single fluid or a single scalar field with a lagrangian of form C = C(</>, d^cbd^cj)) 
^ . it has been proved [k} (see also [§] and that the dark energy perturbation would be divergent as the equation of 
state (EOS) w approaches to — 1. This "no-go" theorem forbids the dynamical models widely studied in the literature 
with a single scalar field such as quintessence, phantom and k-essence to make the EOS cross over the cosmological 
constant boundary. The quintom scenario of dark energy is designed to understand the nature of dark energy with w 
across — 1. The quintom models of dark energy differ from the quintessence, phantom and k-essence and so on in the 
determination of the cosmological evolution and the fate of the universe. 

To realize a viable quintom scenario of dark energy it needs to introduce extra degree of freedom to the conventional 
theory with a single fluid or a single scalar field. The first model of quintom scenario of dark energy is given in Ref. Q 
with two scalar fields. This model has been studied in detail later on. In the recent years there have been active 
studies on models of quintom like dark energy such as models with high derivative term[l2j], vector field (l3l|. or even 
extended theory of gravity [3| and so on, see e.g. [T^ . 

In this paper, we propose a new type of quintom model inspired by the string theory. We will demonstrate in this 
paper our model can realize the equation of state crossing —1 naturally. This paper is organized as follows. In section 
2 we present our model and study its properties especially on the conditions required for the model parameters when w 
crosses over —1. By solving numerically the model we will study the evolution of the equation of state. Moreover, we 
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study the stability property of this model by investigating the quadratic perturbations. The section 3 is the summary 
of our paper. 



II. OUR MODEL 



Our model is given by the following action 1 



S 



d X\J—g 



-V{4>) y/1 - a'V^V^ + 



(1) 



This model generalizes the usual "Born-Infeld- type" action for the effective description of tachyon dynamics by adding 
a term (/>□</> to the usual V M 0V M in the square root. It is known that the 4D effective action of tachyon dynamics to 
the lowest order in V M (/>V M (/> around the top of the tachyon potential can be obtained by the stringy computations for 
either a D3 brane in bosonic theory fl6L Il7| or a non-BPS D3 brane in supersymmetric theory [if 



To this order, we 

have no need to include the operator </>D0 in the action since it is equivalent to the usual V^V^ term. However, we 
cannot exclude its existence in an action such as the "Born-Infeld-type" one when incorporating an infinite number of 
higher derivative terms since now the two terms are in general different dynamically. For example, we cannot simply 
replace the tpUcj) term in the action (TTJ) by the V ' ^V^cf). Further this term in the above generalized action has new 
cosmological consequence as will be shown in this paper. Including an infinite number of higher derivative terms that 
has a significant cosmological consequence has also been discussed in the context of p-adic string recently in (l9j . The 
two parameters a' and (3' in (fT]) can also be made arbitrary when the background flux is turned on [20 ]. Without 
further reasoning, we will take the action (TTJ) as our starting point of phenomenological study. 

As will be demonstrated, the (3' term in (JTJ) is crucial to realize the w across —1. There we have defined a' = a/M 4 
and j3' — (3 /M with a and (3 being the dimensionless parameters respectively, and M an energy scale used to make the 
"kinetic energy terms" dimensionless. V(cj)) is the potential of scalar field (j> (e.g., a tachyon) with dimension of [mass] 4 
with an expected tachyon potential behavior in general, i.e., bounded and reaching its minimum asymptotically, unless 
specifically stated. Note that, □ = -^=0^^/— gg^ LV d u , therefore, in (JXJ) the terms V^V 1 ^ and 0D</> both involve two 
fields and two derivatives. 

The model with operator tpDcj) for the realization of w crossing —1 has been proposed in [121 ] . However in general for 
a model with lagrangian as a sum of operators with a polynomial function of the scalar field <fi and its derivatives, the 
operator </>□(/> can be rewritten as a total derivative term which makes no contribution after integration and a term 
which renormalizes the canonical kinetic term as discussed above. So if one considers a renormalizable lagrangian, 
the operator cf>U(j) will not be included. Ref.[l2| considered a dimension-6 operator as {Ucf) 2 . However in the present 
model, the operator (pUcj) appears at the same order as the operator V^VM) does in the "Born-Infeld-type" action. 
As discussed above, this model appears more natural than the one used in [12 ]. 

With (TTJ) we obtain the equation of motion of the scalar field 0, 



(3 ,V<t>, „ TVo 



2/ 



□ = 



(2) 



where / = ^/l — a'V^V^ + /3'0D0 and = dV/d(j). Following the convention in Ref. [121 ] , the energy- momentum 
tensor is given by the standard definition: 5g ti „S = — J d 4 x^^-T fi,y 5g l j, v , 



(3) 



where -0 = 



W 

dn<f, ~ 2M 4 /' 

For a flat Friedmann-Robertson- Walker (FRW) universe and a homogenous scalar field <fi, the equation of motion 
in |2]) can be solved equivalcntly by the following two equations 



(3V 2 (p M 4 



\M 4 i) 2 j3(j> 
M 4 V> 



/3 2 (f> 2 4M 4 ip 



2a ■ ■ (3V4> 



(4) 
(5) 



1 We have adopted signature (+, — , — , — ) in this paper. 
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where we have made use of the ip as defined before and the first equation above is just the defining equation for ip in 
terms of and its derivatives. H — a /a is the Hubble parameter. 

One can see from equations above the (3 term plays a role in determining the evolution of the scalar field 0. We 
can read the energy density from © as 



and similarly the pressure 



P = Vf + ^W*) + £?ffi-2te, (6) 



With the Friedman equations H 2 = ^j^-p and H = — 47rG(p + p), we now study the cosmological evolution of 
equation of state for the present model. Given w = p/p, we have p + p = (1 + w)p. To explore the possibility of the 
to across —1, we need to check if 4z(p + p) 7^ can be held when to — > — 1. 

Using ([6]) and ([7]) as well as making use of the defining equation for ip, we have 

aV(<f>)<P 2 -d 0V(0) w 

p+p= ^r + ^dt [ ^Fr h (8) 



and from which, 



Eq. © implies that we have either (i)<f> = or (ii)- 2-1220 = /3f P™] when to -> -1. 

Let us assume = first when to — > — 1. With this and Eq. (|9|), we have Jj(p + p) = /30 ^ [ ^ ] = 

— 2^rp/3 2 2 0^D0- Therefore the conditions for having the to across over —1 are (i-l)(f> 7^ 0; (i-2)0 7^ 0; (i-3) JjD0 7^ 

in addition to the = 0. Since the characteristic behavior of V(<p) is to have a finite maximum value at a finite 
and to reach its (exponentially) vanishing minimum asymptotically 2 where one expects = and ^^0 = but 

a non-zero 0, so realizing the above crossing over conditions must happen before reaching the potential minimum 
asymptotically. This implies once the crossing over conditions are met, the field must continue to run away as it 
should be since we have 7^ 0. 

Let us turn to the second case (ii) when to — > — 1. We now have — = pA.(^-^ = j3(f){j-) + /30^(y) which is 
equivalent to saying Y = (a + /3)y0 + /30^(y) = 0. For now, we have ^0. Then the above can be further expressed 
as 2j [(a + /3)ln0 + /31n(y)] = if 7^ 0. 3 Then from ©, the further condition for crossing over is + p) 7^ 
which implies Y ^ or j^[(a + 0) ln0 + /Jlnfj)] ^ °- In 

summary, the following conditions are required for the 
crossing over: ii-l)Y = and ii-2) Y 7^ in addition to the 7^ and 7^ 0. Given the definition of the above Y and 
the characteristic behavior of V(<fi) discussed previously, one expects that both Y and Y vanish asymptotically while 
can reach a non-vanishing value. This implies that the crossing over if occurring at all must occur before the V(0) 
reaches its minimum asymptotically as anticipated. 

Before we demonstrate numerically that the crossing over of the present model can indeed be realized using specific 
examples, let us remark one salient feature of the present model that the (3 term in the present model JT]) is the key 
for realizing the crossing over. One can check from (JSJ that if p = the uj — * — 1 is possible only for = 0. Then 
from ([9]), we will have ^(p + p) = 0, the impossibility of crossing over. 

The analysis above shows various possibilities of our model in realizing the EOS to crossing —1. Now we consider 
some specific examples for numerical calculations of the evolution of the EOS. In Figure [T] we take V{<j>) — Voe~ x ^ 
motivated by string theory and plot the behavior of EOS. In the numerical calculations we have normalized the values 
of the scalar field and Vq, respectively, by the energy scale M. In Figure [1] our model predicts the EOS crossing — 1 
during the evolution and a big-rip singularity for the fate of the universe. Numerically we have checked that = 
when to crosses over the cosmological constant boundary. 



2 This is the behavior of a tachyon potential. 

3 We can always choose <f> 7^ and the = will give a degenerate case which will not be considered here. 
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In Figure [H we take a different potential for numerical calculations. One can see that the EOS crosses over — 1 
during the evolution. When we take (3 positive Figure 3 shows the EOS starts with w < — 1, crosses over —1 into the 
region of w > —1, then transits again to w < — 1. 




FIG. 2: Plot of the evolution of the EOS as a function of In a. In the numerical calculation we take V(<j>) = eA ^^"_ A ^, , and 
Vo — 0.5. For the model parameters we choose A = 1, a = 1, and f3 = —0.8. For the initial conditions we take 4>i — 0.9, 
0i = O.6, (□0) i = £(Cty) i = O. 



Now we study the stability property of this model by considering the quadratic perturbations. Consider a small 
perturbation tt around the background <f>, 

0->#t)+7r(*,x) , (10) 

where the background field <fi in the FRW cosmology is spatially homogenous. Working with the action in JTJ) after a 
shift of the field in (fit))) and a tedious calculation we obtain the terms for the quadratic perturbations 

+...]. " " (ii) 

Interestingly we notice that due to the positivity of the term f~ 3 4> 2 in this model if the coefficient of (V7r) 2 
is positive, the term in front of 7r 2 guarantee to be positively valued. The sound speed characterizing the stability 
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FIG. 3: Plot of the evolution of the EOS as a function of lna. In the numerical calculation we take V((f>) = e \<p^°-\<p , and 
Vo — 0.5. For the model parameters we choose A = 1, a = 1, and (5 — 0.8. For the initial conditions we take 4>i = 0.9, (f>i = 0.6, 

property of the perturbations is given by 

If the c 2 lies in the range between and 1 and the coefficient of tt 2 keeps positive, our model will be stable. In Figure 
2] and Figure [5] we plot the c 2 and the coefficients of tt 2 respectively for the models above in the figures 1-3. From 
these figures one can see for these models c 2 are in the range of 0-1 and the coefficients of 7T 2 are positive, so our 
models are stable. 



III. CONCLUSION AND DISCUSSION 

The current cosmological observations indicate the possibility that the acceleration of the universe is driven by dark 
energy with EOS across —1, which if confirmed further in the future will challenge the theoretical model building 
of the dark energy. In this paper we have proposed a string-inspired model of dark energy through modifying the 
usual effective "Born- Infeld- type" description of tachyon dynamics. As shown in the present work, this modification 
by including a (3 term in the action ([1]) is the key for the EOS crossing —1 during the evolution. 4 Compared to other 
models with w across — 1 in the literature so far the present one is also economical in the sense that it involves a single 
scalar field such as a tachyon and has a motivation inspired from string theory consideration since the new features 
of this model for the dark energy could also be present in the realistic string theory if the whole tower of the higher 
derivative terms is fully included in the effective action. 
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FIG. 4: Plots of the sound speeds (from left to right) for the three models considered in this paper for the numerical calculations 
shown in Figure 1-3. 




FIG. 5: Plots of the coefficient of 7T 2 (from left to right) for the three models considered in this paper for the numerical 
calculations shown in Figure 1-3. 
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